Efficiency of a stirred chemical reaction in a closed vessel 
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We perform a numerical study of the reaction efficiency in a 
closed vessel. Starting with a little spot of product, we compute the 
time needed to complete the reaction in the container following an 
advection-reaction-diffusion process. Inside the vessel it is present 
a cellular velocity field that transports the reactants. If the size of 
the container is not very large compared with the typical length of 
the velocity field one has a plateau of the reaction time as a function 
of the strength of the velocity field, U. This plateau appears both in 
the stationary and in the time-dependent flow. A comparison of the 
results for the finite system with the infinite case (for which the front 
speed, Vf, gives a simple estimate of the reacting time) shows the 
dramatic effect of the finite size. 
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Numerous physical, biological and chemical systems show 
the propagation of a stable phase into an unstable one [|l]j2|]. 
When this phenomenon takes place in a fluid, one generally 
speaks of front propagation in advection-reaction-diffusion 
(ARD) systems. Under this generic name one indicates many 
different processes, e.g., the propagation of plankton popula- 
tions in ocean currents [|^], the transport of reacting pollutants 
in the atmosphere (e.g. ozone) [Q], or the premixed combus- 
tion [|]. 

In the last years much effort has been done to study the 
influence of an advection field on the front dynamics. In par- 
ticular, it is well established that the front speed in a laminar 
or turbulent fluid is enhanced with respect to the propagation 
in a medium at rest In the context of (premixed) com- 

bustion processes the flame front area is proportional to the 
front speed and, therefore, an increasing of the front area due 
to the fluid stirring gives rise to an enhancement of the burning 
efficiency, that is, the system burns faster. 

It is important to note that most of the theoretical studies 
and, in particular, the above results about the enhancement of 
burning efficiency, have been shown for an infinite-size sys- 
tem (in the propagation direction). Moreover, in order to in- 
troduce well-defined mathematical quantities one is forced to 
work with an infinite (or with periodic boundary conditions) 
system. This is the case of the front speed, which is an asymp- 
totic quantity well defined only for an infinite system. 

However, from a practical point of view one usually has to 
treat cases where the size of the domain is not much larger 
than the typical length scale of the velocity field [^j]. The 
spreading of organisms in a lake or in a small closed sea basin, 
or the combustion of fuel in a machine motor are two clear ex- 
amples where this may happen and, therefore, non-asymptotic 
properties can be very important [P,p^. 



In this work we treat the case of an ARD process confined 
in a closed area. Beginning with a small quantity of material 
in the stable phase (in the following called burnt material), we 
numerically compute the time needed for a given percentage 
of the total area to be also burnt (called in the following, the 
reacting or burning time). The velocity field is of cellular flow 
type, that is, formed by circulating cells of fluid. Both, the 
stationary and the chaotic time-dependent cellular flow will 
be considered. Our main result, obtained either for the time- 
independent and for the time-dependent flow, is that increas- 
ing the typical velocity of the field one has a saturation of the 
burning rate. This saturation happens when the advection time 
scale is much faster than the reaction time scale. Also, we 
compare our results with the infinite-size case, studying the 
crossover from finite size systems to the asymptotic regime. 
We observe that the relevance of the system size is more im- 
portant than expected a priori. 

Let us consider the simplest non trivial case described by a 
scalar field 0(x, t) which represents the concentration of reac- 
tion products, such that 9 is equal one in the space-time co- 
ordinates where the reaction is over (the stable phase), and 9 
is zero where there is fresh material (the unstable phase). The 
dynamical evolution of this field is given by 



dt9 + v-V9 = DV^e 
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where v(x, <) is the two-dimensional velocity field, D is the 
diffusion coefficient, and f{6) is the reaction term, where r is 
the time scale for the reaction activity. For the reaction term 
we use the Fisher-Kolmogorov-Petrovskii-Piskunov (FKPP) 
nonlinearity JTl]], f{9) = 9{l — 9). Concerning the velocity 
field, V, we first adopt a simple stationary incompressible two- 
dimensional flow defined by the stream function 
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being the parameter U the maximum vertical velocity of the 
flow, and L the size of one cell. For a study of the trans- 
port properties in the field (||) see ref. [p^; the asymptotic 
behaviour of front propagation is discussed in ref. [pj|]. The 
equations of motion for a fluid element are given by 
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In this work the reaction processes described by (|l]) take place 
in a closed recipient. This confinement is implemented by 
assuming rigid boundary conditions on the box < y < L 
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and < a; < nL, where n is the number of circulating cells 
of the flow. One approaches to the asymptotic case increasing 
the value of n. 

The settling of the problem is completed when we indi- 
cate the initial conditions, i.e., an initial spot of burnt material 
which starts the reaction. Thus, we use in all our numerical ex- 
periments a small circle of radius r filled with stable material 
{9 — 1), that is placed at the initial time in the box filled with 
unstable material {6 ~ 0). The initial coordinates of the center 
of this circle are {x = r,y = L/2) (the circle is on the border 
of the box; this mimic the injection of reacting material from 
the outside). As anticipated, the principal observable under in- 
vestigation is the time needed for a given percentage of the to- 
tal area to be burnt. We define S{t) = J^^ dxdy9{x, y, t) 
as the percentage of area burnt at time t, where AS — L^n 
is the total area of the container. In our case, by choosing an 
appropriate r the initial burnt material is S{Q) — O.OOS/ti, 
which is the 0.5% of one cell. The reacting or burning time ta 
is defined as the time needed for the percentage a of the total 
area of the recipient to be burnt, i.e., S{ta) — a. 

Numerically, to integrate (|l]) we use the Feynman-Kac (FK) 
or stochastic Lagrangian approach. In this algorithm the field 
evolution is computed using the Lagrangian propagator plus 
a Montecarlo integration for the diffusive term. Then, the re- 
action propagator accounts for the reacting term (for details 
see [ pj[p^ ). We also impose a rigid wall condition in the 
boundaries, in order to avoid that any fluid particle leaves the 
container, which could happen due to the noise term added to 
the velocity field in the Lagrangian approach. 

We first show in Fig. (|l]) the influence of the velocity on the 
reacting efficiency, when different percentages, a, of the final 
burnt area are considered. 




FIG. 1. ta against the flow strength U with r — 0.4, n — 1 and 
for various a: a = 0.2 (+), a = 0.5 (x), a = 0.7 (□), a = 0.9 
(o). 

Increasing the velocity of the flow, U, decreases mono- 
tonically until a plateau is reached. Then, a further increasing 
of the flow velocity (U > 15, similar for different a's) does 
not decrease the burning time We remark that this effect 
also appears for different finite values of the system size (dif- 
ferent n's), and different chemical rates r. At first glance, the 



appearance of the plateau seems to be surprising: in a closed 
container and for very high stirring intensity, increasing fur- 
ther more this intensity there is not an enhancement in the 
burning front propagation, that is, one does not improve the 
efficiency of burning. 

The existence of the plateau can be understood noting that 
it is reached only when the reaction time, r, is large com- 
pared with the advection time, Tq = L/U . In this case, in 
the first stages of the process the active material invades the 
whole container because of advection and diffusion. Then the 
reaction term is the final responsible for the cell burning. 

A direct comparison between the finite and the infinite sys- 
tem is interesting. In Fig. (||) we show the burning time scaled 
with the system size, i.e., ta/n, against the typical flow ve- 
locity, U, for some values of n. We also we plot the data 
obtained for an infinite system, which have been calculated 
from the front speed f / of the infinite system data according 
to 

nL 

^ — , (4) 
which is expected to hold for a close to 1 and large n. 
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FIG. 2. The burning time per unit-cell, ta/n, at various U for 
r = 0.4 and a — 0.9. The plots are for different system's sizes: 
n = 1 (+), n = 2 (x), n = 4 (*), n = 8 (□) and n = 12 (o). It 
is also shown (•) the burning time calculated using the front speed: 

t = nL/vf. 

Figure (||) shows that the asymptotic reacting time (given 
by Vf) is reached only in the large size limit, i.e., n large, 
while the dynamics of small systems is dominated by the non 
asymptotic properties of the evolution. 

Since the considered bidimensional velocity field (^ is sta- 
tionary the Lagrangian trajectories are not chaotic. Neverthe- 
less also in the case of Lagrangian chaotic trajectories, ob- 
tained with a time-dependent flow, the burning time shows the 
same qualitative behaviours shown in figures and (|). Let 
us consider a time dependent flow whose streamfunction is 



ip{x, y,t) — U sm{x + B[x) cos{ujt)) sin{y) . 



(5) 



This is sufficient to induce Lagrangian chaos [ |16[ | in the evolu- 
tion of passive tracers advected by the velocity field (^ gen- 
erated by (Bb. Because we are dealing with closed systems 
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B{x) is constructed in such a way that it is zero near the 
boundary of the system and almost constant, Bq, otherwise: 

B{x) = Bo{l - exp(-l/a;) - exp(-l/(nL - x))). 




FIG. 3. ta against the flow strength U for various a in the un- 
steady case (B = 1.1 and u = 2.09): a = 0.2 (+), a = 0.5 (x), 
a = 0.7 (□), a — 0.9 (o). The flow is confined in two cells. 

In Fig. we show the curve ta against U for different val- 
ues of a. At difference from the previous case, when unsteady 
flow is concerned there is not a simple plateau in t^, but an 
oscillatory behaviour due to the interplay between the oscilla- 
tion period of the separatrices and the circulation time inside 
the cell. It happens that circulation and oscillation "synchro- 
nize" producing a very efficient and coherent way of transfer- 
ring passive particles from one cell to the other A similar, but 
much more impressive, feature occurs for the effective diffu- 
sion coefficient in the horizontal direction [|l7|] . Anyway, for 
low values of U the mixing induced by the time dependence 
makes that the system burns quite faster. For higher U, the 
mixing properties of the flow are not sufficient to improve fur- 
thermore the burning efficiency, and at the end, the typical 
time for the cell burning is proportional to the reaction time- 
scale. 

Summing up, the physical mechanism of reaction in 
a closed recipient for both time-independent and time- 
dependent cellular flows are quite similar When the typical 
time-scale of the velocity field, L/U, is larger than the reac- 
tion time-scale, r, the initial condition is quickly spread along 
the whole system due to advection (and diffusion), in such a 
way that in all the container the value of 9 is small but differ- 
ent from zero. Then, basically only reaction plays an impor- 
tant role, and a further increasing of the flow velocity does not 
decrease the burning time. 

Let us study the dependence of the saturation time ^^(t), 
which is the value of the burning time in the plateau, as a 
function of r. In the unsteady case, we choose the minimum 
value of ta at varying U as the saturation time. 

In Figure (Q) we show the result for the unsteady case, 
which can be interpreted following the same arguments used 
to explain the existence of the plateau. Essentially the burning 
process in the case of large U can be divided in two steps, the 
initial mixing regime and the reacting dominated regime. In 



fact, from the initial condition (in which only a small fraction 
of the volume is active) one has a rapid spreading (say in a 
time t) because of advection and diffusion. After the spread- 
ing one has an exponential increasing of the field due to the re- 
action term. The time i is expected to be a decreasing function 
of U (with a limiting value for sufficiently large U). This im- 
plies, together with a dimensional argument, (r) = t + baT. 
Thus, there is a linear dependence of t^ with r as shown in 
Figure (^, which has been obtained for the time-dependent 
flow, but similar results hold also for the steady case. 
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FIG. 4. Plot of ^^{t) vs t at various percentages of filling a in 
the time-dependent case (B — 1.1 and uj = 2.09): a — 0.2 (+), 
a = 0.5 (x), a = 0.7 (□), a = 0.9 (o). For each curve we have 
superimposed the linear fit t„(r) = t + baT, where t and 6^ are 
fitting parameters. Here n = 2. 

The value of the slope of the saturation time against t, ba, 
can be analytically studied. We have that, in the regime of 
very high U such that the plateau is reached, after the time 
t the initial condition is spread out through the whole con- 
tainer and one can approximate 9{x, y, t) ^ d{t), being 9{t) 
a rough average of the 9 field in the container, which evolves 
following only the reaction part of (|l]). This is because, as we 
have argued before, after t the important physical mechanism 
of burning comes from the chemical activity and not from the 
mixing due to the advection and diffusion. Then S{t) = 9{t) 
for t > i, and so: 



at T 



(6) 



This can be integrated from i to t, taking into account that we 
can approximate S{t) ~ 5'(0), i.e., the initial condition is just 
spread out in the system in the initial stages of the process. 
One has that 



log 



S{t) \ _ / ^(0) \ ^ t_i 
l-Sit)J Hi -5(0)7 r 

Finally, as S{t^) — a we get 

, ~ , /a(l-S'(0))\ ~ , 



(7) 



(8) 
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which gives the dependence of ba on the percentage of burnt 
material, a. 

For the time-independent cellular flow the reasonings fol- 
low closely the former ones. In fact, in the regime of large 
U, when all over the cell (by diffusion) there is also a small 
quantity of active reagent, the reaction process can begin with 
an averaged S{t) (in a mean-field sense). 

Despite the numerous approximations done to obtain 
it is in excellent agreement with the numerical results, see 
Fig. (|]), confirming once again the physical mechanism we 
think that give rise to the existence of the saturation time. 




FIG. 5. The slope, ba, of the saturation time against r, that 
is, the slope of the curves in Fig. (with an additional value for 
a — 0.95), vs the percentage a. With (□) the numerical values 
for time dependent flow, with (*) the numerical value (rescaled) for 
steady flow, and with the solid line the prediction given by (^. 

Summarizing, we have performed a numerical study of an 
advection-reaction-diffusion system confined in a closed ves- 
sel, using stationary and time dependent cellular flows. Be- 
ginning with a small quantity of the active phase, we have 
calculated the time needed for a percentage of the total area 
to be burnt. Thus, our numerical experiments may represent 
the spreading of an organism in a lake or the combustion of a 
material in a vessel. The main lesson to learn from our stud- 
ies is that the influence of the system size has been shown 
to be very important JT^ . In particular, we have shown that 
for very large flow velocities the reacting time saturates, giv- 
ing rise to the unexpected result that increasing further more 
the flow velocity there is not a decreasing of the time to burn 
the material. We have to mention that a similar scenario, i.e., 
the appearance of the plateau in the burning time, has been 
obtained for other types of chemical reactions f{6), like the 
Arrhenius f{9) = (1 — 9) eiq){-0Q/9) (9o constant) or the 
Zeldovich function f{9) = 9"'{1- 9), with m = 2. 
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